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Quantum walks subject to decoherence generically suffer the loss of their genuine quantum feature, 
a quadratically faster spreading compared to classical random walks. This intuitive statement has 
been verified analytically for certain models and is also supported by numerical studies of a variety 
of examples. In this paper we analyze the long-time behavior of a particular class of decoherent 
quantum walks, which, to the best of our knowledge, was only studied at the level of numerical 
simulations before. We consider a local coin operation which is randomly and independently chosen 
for each time step and each lattice site and prove that, under rather mild conditions, this leads to 
classical behavior: With the same scaling as needed for a classical diffusion the position distribution 
converges to a Gaussian, which is independent of the initial state. Our method is based on non- 
degenerate perturbation theory and yields an explicit expression for the covariance matrix of the 
asymptotic Gaussian in terms of the randomness parameters. 
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I. INTRODUCTION 



Quantum walks describe the time evolution of a single quantum particle with internal degrees of freedom, for which 
both space and time are discrete. We study here the case where the underlying space is assumed to be an infinite 
lattice of arbitrary dimension. The dynamical rule is given by a unitary operator composed of a coin operator acting 
on the internal degree of freedom only, in a generally site-dependent way, and a fixed shift operator translating the 
particle by finitely many lattice sites depending on its internal degree of freedom. We are interested in a situation 
were the coin is varied randomly, as a way to model the imperfections of experimental realizations. If the distribution 
of coins is highly peaked around a fixed one, we would expect to see a coherent walk with a linear increase of the 
standard deviation with the number of time steps t, at least for some time. In the long run, however, the randomness 
will be felt, and it is this regime we will study. Depending on where we put the random dependence we can distinguish 
four cases, summarized in the following table, and discussed in turn below. 



coin 


spatially 


dependence 


fixed 


random 


temporally 


a ~ t 


a ~ 1 


fixed 






temporally 


a ~ \ft 


a ~ \ft 


random 







Table I. The type of randomness of the coin operator determines the dependence of the standard deviation a on the number 
of time steps t of the quantum walk. The cases where the random coin operator is fixed in time respectively space have been 
studied in the literature and it was shown that they lead to localization (a ~ 1) respectively diffusive behavior (a ~ \ft). Our 
aim is to analyze the case where both types of randomness are combined and we prove that the asymptotic behavior of such 
quantum walks is diffusive. 



Coherent walks, i.e., walks without any randomness, have been found to be useful in search-like algorithms |Kem05[ 
IAmb03|. IAmb07| ICCJY091 IFGG08] . precisely because they spread faster than classical walks [Amb03j . which have a 
similar algorithmic use. They are also the simplest case of quantum simulators, since they can be understood as the 
one-particle sector of so-called quantum cellular automata JSWj , which are quantum systems on a lattice of infinitely 
many interacting quantum particles. There has also been don e a lot of experimental work to implement quantum 
walk s in a var i ety of physical setup s, starting wi th cold at oms |KFC + 09] and followed by experiments with trapped 
ions [SMS + 09UMSE+llllZKG+10] and photons [SCP+10] , 

Temporal fluctuations are implemented by a walk operator that is random in time, which means that for every time 
step a different walk operator has to be applied, keeping, however, the spatial translation invariance in each step. It is 
clear that we have to take the expectation value over all possible sequences of time dependent coin operators in order 
to model fluctuations of the coin parameters. This expectation value turns the formerly unitary time evolution into 
a decoherent quantum channel. This model has been studied in a number of examples [CSB07, ADSS07, SBBH03, 
IBCA021 ISK081 IKBH06j and in great detail in |AVWW111 |Joyll| IHTTT] and it was shown that such a time evolution 
generically leads to diffusive behavior of the quantum walk, which means that the standard deviation of the position 
probability distribution grows proportionally to the square root of the number of time steps. 

Spatial fluctuations of experimental parameters correspond to the case where the time evolution is still unitary, 
and the same unitary in every step, but the coin operator is random in space. For continuous time, i.e., Hamiltonian 
systems this is the well known Anderson model of disordered crystals, which exhibits localization. This means that 
the Hamiltonian almost surely has purely discrete spectrum, and the position distribution does not spread at all. 
The case of quantum walks on a one-dimensional lattice subject to spatial disorder has been studied in a number 
of examples both numerically |LKBK101 IQKllj and theoretically |Kon09bl IKon09al IJM101 lASWllj . These results 
show that, at least in one-dimensional systems, spatial disorder implies dynamical localization, meaning that after 
arbitrarily many time steps the quantum walker is confined to a finite region of the lattice, up to exponentially small 
corrections. 

In this paper we examine the case wher e both typ es of disorder appear simultaneously. Such quantum walks have 
been studied numerically for example in |RSA + 05] and [LKBK10] and the simulations indicate diffusive behavior. 
We use the coin and shift decomposition mainly to have a precise meaning for independently identically distributed 
randomness. For this setting we develop a general theory of the asymptotic position distribution and find diffusive 
scaling. More precisely, the scaled limiting distribution is exactly gaussian and independent of the initial state. This 
distinguishes the present case from only temporal randomness, where we get gaussianness only in every momentum 
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component. Since momentum is conserved, a residual dependence on the initial state remains, and since the diffusion 
constant depends on momentum, the resulting mixture of Gaussians is no longer a Gaussian. 

Let us briefly outline the structure of this paper. We start in section [TT] by the mathematical formulation of the 
model, followed by the general examination in section III We continue our discussion by the application to a variety 
of examples in section |IV| In section [V] we comment on generalizations of our results to the case of more general 
quantum walks and we conclude in section |VI| by discussing some open problems left for future research. 



II. MODELS OF QUANTUM WALKS WITH RANDOM COIN 

Abstractly, quantum walks can be defined as a discrete time evolution of a quantum particle with internal state 
space JC moving with strictly finite propagation speed on a lattice Z s . Usually, one also assumes translation invariance 
of the time evolution which then yields a structure theorem [AVWWTT] for the class of all translation invariant and 
possibly decoherent quantum walks. Hence, the underlying Hilbcrt space is H — I 2 {I s ) £g) fC and quite commonly a 
single time step is realized by a composition of a local coin operator C and a conditional translation operator called 
(state-dependent) shift operator S. Throughout this paper we will assume that the shift operator S is given by a 
unitarily implemented quantum channel, hence, if C is also unitarily implemented we can represent the quantum walk 
by a unitary matrix W acting on £ 2 (Z S ) ® JC. We denote the unitary operators corresponding to the coin respectively 
shift again by C respectively S, and hence, W — S ■ C . To begin with we specify S and C in the case where both are 
unitarily implemented. 

We denote the elements of the standard basis of £ 2 (Z S ) <Ei JC by \x <E> i), where x E Z s labels the positions and 
i = 1, . . . , dim JC labels a basis of JC such that S is given by 

S\x®i) = \x + Vi®i), (1) 

with some vectors Vi € Z s . A single time step at time t is generated by the unitary operator 

W t = S'C t , c t = U Xtt , 

where U x t is a unitary matrix of dimension dim/C depending on the time t and the lattice site x. Ideally, the coin 
operator is translation invariant and constant in time, that is U x ^ — U . In this case we can write C — <£> U and 
it is well-known that this generically leads to ballistic behavior of the quantum walk, that is, the standard deviation 
of the position distribution grows linearly with the number of time steps a(t) ~ t. 

If the coin operator at a fixed time t is translation invariant, but varies in time, we have U Xx t = Ut- In this case, 
the interpretation of fluctuating coin parameters corresponds to a lack of controllability of the unitary Ut- In other 
words, instead of a deterministic sequence U\, . . . , U t applied sequentially to an initial state we actually have to take 
the expectation value over all possible sequences of time dependent coin operators. In fact, we cannot control which 
coin operator happens at a certain time and according to quantum mechanics and its statistical nature we have to 
repeat the experiment several times, each with a different sequence of coin operators, in order to extract the position 
distribution of the quantum walk after a fixed number of time steps. Let us assume that the coins U t are distributed 
independently and identically in time according to some measure v on U(JC), the space of unitary operators on JC. 
We identify the underlying probability space £1 with U (JC) and an element u € f2 uniquely determines an operator 
U u £ £/(7C) |Not] . With this notation we can describe a single time step of an observable A in the Heisenberg picture 
by the application of a decoherent quantum channel according to 

W(A) = S* • E [(lz. ® U*) ■ A ■ (1 Z3 ® UJ)] ■ S 

UJ 

where the expectation value is taken with respect to the probability distribution v of the coin operators 6 U(K). 
This case has been studied in |AVWWlT| |Joyll[ IHJ11] and it was shown that such a time evolution generically leads 
to diffusive behavior, by which we mean that the standard deviation of the position distribution grows proportional 
to the square root of the number of time steps o~(t) ~ \ft. 

If on the other hand the coin operator is constant in time, at least for a large number of time steps, but inhomoge- 
neous in space, we have U x j — U x . Here, the time evolution after t time steps is generated by the t-fold concatenation 
of W = SC, where C = 0^ U x = J2 X \ x )( x \ ® U x is the coin configuration generated by the unitary matrices U x . If 
now the fluctuation of the experimental parameters is only spatial one gets a unitary time evolution, but if on the 
other hand there is a temporal fluctuation in C on a large time scale, comparable to the duration of a single run of 
statistical data collection, one needs to take the expectation value over all possible spatial realizations of the coin 
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operator C . Assuming again that the matrices U x are distributed independent and identically according to some 
measure v on U(IC) we can write the time evolution after t time steps as 

W t (A)=M 

u 

where the expectation is taken with respect to the spatial configurations of C, which mathematically corresponds 
to the infinite product measure v°° of v defined on the probability space W^IA^IC). The one-dimensional case 
H = i 2 {%) <g> C 2 was analyzed in |ASW111 IJMlOj and it was shown that such a time evolution yields Anderson 
localization, that is, up to exponentially small corrections, the position distribution of the quantum walker has finite 
support on the lattice Z under rather general assumptions on the distribution v. Those results also apply to the case 
where we do not average over all possible spatial configurations, in other words, almost all possible configurations 
already show Anderson localization and large-scale temporal fluctuations are not required for localization. 

In this paper we will analyze the combination of the former mentioned cases. Fluctuations of the coin parameters are 
now assumed to happen in space as well as in time on the scale of a single time step. Mathematically, this means that 
the i-fold time evolution of a single run of the experiment is given by the unitary operator W t ■ ■ ■ ■■ W\ — SC t ■ ■ ■ ■■ SC\ . 
Similarly to the other models we have to take the expectation with respect to the distribution of the coin operators, 
but now in space and time. A crucial assumption we impose on our model is again that the coins are independent 
and identically distributed both in time and space according to a measure v on U(IC). Consequently, one time step of 
the evolution can be written as 



S* ■ @ C/* (x) • A ■ 




W{A) = S*C(A)S, (2) 

where S denotes the unitary shift operator and C the decoherent coin operator stemming from the expectation value 
with respect to all possible coin realizations. The action of the averaged coin operator C on a generic operator 
A = J2 X y \ x )(y\ ® A xy can be written as 

C Ij2\x)(y\ <8> A xy I = J2 \ x )(v\ ® Uxy J y{duj)U* u A xy U u (3) 

\x,y / x,y ^ 

+ (1 - 5 xy )U*A xy uS) , U = ( v(du)U u . 



Now since the distribution v is independent of the lattice site x it follows that W itself is a translation invariant 
operator. This is similar to the model considered in |AVWWTT] . the crucial difference being that there the existence 
of a Kraus decomposition of W in terms of translation invariant Kraus operators was assumed. This, however, is not 
the case in this model, where a Kraus decomposition is given by the Kraus operators corresponding to all possible 
realizations of coin operators, which is a decomposition into non-translation invariant Kraus operators. We will further 
develop our method used in |AVWWTT] in order to cope also with the case of fluctuations in space and time and 
prove that this in fact leads to diffusive behavior. 



III. THE PERTURBATION METHOD 



Our goal is to determine the scaling of the standard deviation a(t) of the position probability in time, in particular 
to distinguish between ballistic (cr(t) ~ t) and diffusive (cr(t) ~ behavior. Using perturbation theory we compute 
the asymptotic limit t — > oo of the position distribution. 

We start this section with a description of the general theory we are going to apply to quantum walks according 
to ^ and (|3b. Since our method is based on perturbation theory of infinite dimensional operators we need to verify 
analyticity of the series expansions explicitly. This is done in the second part of this section. The actual analysis of 
our quantum walk model is postponed to the third part of this section, where we summarize our results. 



A. General Theory 



Similarly to the approach in [AVWWTT] we compute the characteristic function of the position distribution of the 
quantum walk with initial state po after t time steps scaled by e 



C^(A)=tr( Po W t (e ieA - ( 3)) 



(4) 
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and determine the limit C(A) = lim Ct e(A), where e is either chosen as X/t in ballistic scaling or \j\ft in dilfusive 

t— >oo ' 

scaling. The method we are going to incorporate is based on perturbation theory of bounded operators Kat95] . The 
idea is to introduce a similarity transform 

<f> : B{H) — > B{%) (5) 
A i— > <f){A) = Ae ieAg 

and define the operator 

W £ := 4>- 1 o W o (6) 
on i.e. W e (A) := W(Ae i£AQ )e- ieAC? , and rewrite 

W 4 (e i£A - Q ) = W^(l)e ieA - Q . (7) 
When inserting this into C(A) we can neglect the factor e leX '® since po is a trace-class operator and hence 

C(A) = lim tr (p W*(l)e i£A Q ) = tr f lim (p Wt(l)e ieXQ )) (8) 

= trf lim (W*(l)) lim (p a e iEXQ ) 

\t— >oo — 1 

= trfp lim W*(l) 



In ([8| we assumed that the limit of W'(l) for t oo exists in operator norm with appropriate scaling of e. In 
fact, our goal is to determine this limiting operator via perturbation theory and by inserting it into ([8]) we get the 
characteristic function of the asymptotic distribution of W and po- We interpret e as a perturbation parameter, so 
we act with a high power of the perturbed operator W e on the eigenvector 1 of the unperturbed operator W = Wo. 

Before we deepen our analysis let us sketch the results to be expected. The operator W e approaches W as e — » 0. 
Moreover, the perturbed eigenvector A e , obeying W e (A e ) — p, e A e , approaches 1 as e — > 0, hence, we expect 

W^(l)-+/4l. (9) 

t— > oo 

Consequently, the characteristic function is given by C(A) = lim^oo p} e , which, in contrast to the case considered in 
[AVWWTT] . is always independent of the initial state po- From the perturbation expansion p e = 1 + fi'e + p"e 2 /2 + . . . 
with p! = iv ■ A and p" = — A T ■ D ■ A we get in ballistic scaling e = 1/t the characteristic function 

C7 1/t (A) = «>\ (10) 

which is the characteristic function of a point mass at v corresponding to a constant drift with velocity v. If p! = 
we can look at the diffusive scaling e = 1/ \ft of W which yields 

o 1/Vi (x) = e - AT D ' A (11) 

corresponding to a position distribution which is a Gaussian with covariance matrix D. 



B. Analytic Perturbation Theory for Quantum Walks 



To begin with let us put this problem rigorously into the context of perturbation theory. It is convenient to consider 
vectors ip — J2 X \ x )( x \ ® V'x € £ 2 {Z S ) ®K as functions tp '■ x ^ m other words, we identify £ 2 (Z S ) ®JC with the set 
^ 2 (Z S ,/C) of all /C-valued square summable functions on I s . Then, a translation by a vector y € Z s on £ 2 (Z S ) ® /C, 
which we denote by T y , can be defined via T y ip : x i-» ipx-y With the help of these T y we define translations of 
bounded operators A € Z?(^ 2 (Z S ) ® /C) by Tj,(A) = T y AT* and denote the set of all translation invariant bounded 
operators on £ 2 (Z S ) Cg) /C by 7ic,s C S(£ 2 (Z S ) ® /C). The defining equation for 7jc, s is 

AeTc s ^ ^(,4) = A VyeZ s . (12) 

Now, let us argue why 7jc s constitutes a vector space on which W e acts, which is expressed by r y o W e = W E o t v 
for all y G Z s . Indeed, though the similarity transform on 6(^ 2 (Z S ) eg) /C) does not preserve translation invariance, 
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more precisely r y o <f> = e lsX " v (f> o t v , the operator W e = o W o <f> commutes with translations T y due to the fact 
that W preserves translation invariance and the two appearing phase factors cancel each other. 

The analysis of translation invariant operators A 6 T/c, s and maps acting on Tic s is much simplified by introducing 
the Fourier transform on £ 2 (Z S ) K. via 

(^)(p) = (27T)"* eiX ' P ^> PS [0,2tt) s . (13) 

Then each A G Tic.s becomes a multiplication operator in Fourier space, i.e. it exists a unique matrix valued function 
p i ^ j4(p) € 2?(/C) such that (FAij))(p) — A(p)ip(p) |SW71| . and vice versa. This leads to the following sesqui-linear 
form on Tk.s 

{m '-=J^y I dSp &Lic trKA * {p)B{p) ' (14) 

[0,2ir) s 



which in fact is even a scalar product turning Tic s m t° a separable Hilbert space. Another way of interpreting (14 1 is 
to observe that a translation invariant bounded operator A G Tk,s is fully characterized by its matrix entries at the 
origin, i.e. if we expand A in position basis as A — J2 X yez s \ x )(y\ ® A xy with A xy g B{K) then = Ao jX _ y . The 
corresponding multiplication operator in Fourier space is now given by A{p) — X^ez^ eip x A x o and therefore we get 
the alternative expression 

( A \ B ) = mhc E ^ A toB x o , (15) 

x£Z s 

which is finite since A and B are bounded operators. We denote the norm on Tic s induced by this scalar product by 
|| . ||, thus, \\A\\ 2 = (A\A) for all AeT K , s . 

Now, after we have introduced the Hilbert space Tc,s we consider the specific form of quantum walks according 
to ([lj, ([2]) and |3j. The shift operator S is represented in momentum space by conjugation with the p-dependent 
dim /C-dimensional matrix 



S(p) 



/ e iVlP ... \ 

'•. 

\ : e i«dimiC-p J 

hence, the operator W acts on a translation invariant bounded operator A(p) in the following way. 



(16) 



W(A){p) - S(p)* U v (du)UZAoU u (17) 
+U*(A(p)-A )u)s( P ) 



In this equation Aq denotes the p-independent term in A(p), which can be represented as Aq = (2ir) s J d s pA(p). 
The modified operator W £ is now given by 



W e (A)(p) = S(p)* l^j v(duj)u:AoU^ (18) 

+U*{A{p)-A )u) S( P + e\). 

where the momentum shift of eX arises from the definition of the perturbed walk operator W e and the fact that the 
coin operator C commutes with e leA < ^, together with the equation 

e l£XQ S(p)e~ leXQ = S(p + eA) . (19) 

Our goal is to apply non-degenerate perturbation theory to the operator W £ , in particular, the aim is to determine 
the first and second order equations of the perturbation theory in e. The correctness of the results obtained by 
equating coefficients of powers of e is in fact non-trivial since W e is defined on an infinite dimensional Hilbert space 
Tk.s- Hence, we first have to establish analyticity of W £ , the eigenvector A e and the corresponding eigenvalue /i e 
with A e — > 1 and /i e —> 1 as s — !• 0. This can be done by using the following theorem, which is an adaption of the 
well-known theorem of Kato and Rellich Kat95, RS78 to our setting. 
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Theorem III.l (Kato-Rellich) . Assume that the operator W £ is bounded on 7jc, s an d the limit of the difference 
quotient 

lim^^ 



A^O A 

which we then call the derivative of W £ at e, exists in operator norm for all e in an open subset of C containing the 
origin. If the eigenvalue equation W(l) = 1 of the unperturbed operator W = Wo is non-degenerate, then, for small 
enough e, there exists an analytic eigenvector A e with analytic and non-degenerate eigenvalue fi E such that A £ — > 1 
and [i e — > 1 as e — > 0. 

For a proof of this theorem we refer to [Kat95, RS78J. Now, we are left to prove that W e is bounded, differentiable 
and the eigenvalue 1 of W is non-degenerate. In order to ensure the non-degeneracy of the perturbation theory we 
assume that some power of W is strictly contractive on the orthogonal complement of 1 . Below we identify a class of 
quantum walks W which are strictly contractive on their own, i.e. ||W(vl)|| < \A\ for all A _L 1. 

The following proposition in conjunction with theorem |M.1| assures the analyticity of the perturbation theory of 
W e . 

Proposition III. 2. Let W be a quantum walk according to equations ^ and ([3| and W e be defined by ([6]). 
Assume that some power of W is strictly contractive on {l} -1 ", that is, there exists n £ N such that ||W™(^4)|| < \\A\\ 
for all A A. 1. Then we have the following conclusions: 

i) The derivative of W e exists in operator norm for all e € C. 

ii) For all s £ C the operator W £ is bounded with ||W e || op < max |e 1£A " Ui | . 

i 

Hi) The eigenvalue 1 o/W is non- degenerate. 
Proof. It is easy to check that the operator W' e defined via its Fourier transform 

W' E (A)(p)=W e (A)(p)-iA, 

where A is the dim /C-dimensional diagonal matrix with matrix entries A • Vi on the diagonal, is indeed the operator 
norm limit of the difference quotient 

r W e+A -W e 
lim 



A^O A 

which proves i). 



By ( 16 ) we have S(p + e\) = S(p) ■ S(e\), which implies 

(W £ (A)\W £ (A)} = (W{A)\W(A)S(e\)S*(e\)) < max |e ieA " u ' \ 2 (W (A)\W (A)) , 

i 

and hence, by the basic definition of the operator norm, we have ||W e ||™ < max \ e leX ' Vi |||W|| op . Now, ii) follows from 

i 

1 1 W 1 1 op = 1, which we prove next. 

By writing W(A) = S*C(A)S and observing that S is unitary we get ||W(A)|| = ||C(A)||, and hence ||W|| op = 
|| C || op . We denote the first part of the coin operator in ^ by T, that is, 



T(A) = J v{dw)U*AU u 

and hence 

C(A)(p) = T(A ) + U*(A(p)-A )U. 

The Hilbert space Tic, a can be decomposed into a direct sum of two orthogonal subspaces To and 7q defined via 

To = {A e Tk, s ■ A = A } (20) 
= {A e T K , S : A = 0} . 

Clearly, C(T ) C To and C(7 J -) C T^, from which it also follows that W(%) ± ^(Tq 1 -). Hence, 

||W|| op = \\C\\ op = max{||T|| op , \\U* .U\\ op } , 



where we consider T respectively U* .U as map on 7o respectively 7~q~- First, let us note that 

\\T(A)\\ 2 < J J v(du)v(du')\tr X :UZA'U u UZ,AU U '\ 

< J J v(duj)v(duj')y/tT K A*Ay/tT K A*A 

in other words, ||T|| op < 1. Now, let A and B be positive operators on K,, then we have 

tncAB< \\A\\ op ■ tT K B , 
and hence, by applying this inequality to tijcC(A)*C(A), we get 

tTicU*A*UU*AU = tT>cUU*A*UU*A (21) 

< \\UU*\\ op ■ tv K A*UU*A 

< \\UU*\\ 2 op -tvjc AA* 
<\\U\\ 4 op -tv^AA*. 

Clearly, ||E/|| op < / v{du)\\U UJ \\ op = 1 which finally proves ||f/*.E/|| ap < 1, and hence, ||C|| op < 1. 

Any eigenvector of W is also an eigenvector of W n with eigenvalue raised to the n-th power, hence, the contractivity 
of W™ yields statement in). □ 

Of course, the contractivity of W or a power of it depends on the probability distribution v of the coin operators U^. 
In particular, the operator W itself is strictly contractive if the operators U^U^, with u;',ui e £1 fulfill the following 
definition. 

Definition III. 3. A set of matrices {Mi : i e I}, where I is an index set and each Mi acts on IC, is said to be 
irreducible if any invariant subspace is trivial, that is, if S is a subspace of K such that M^S C S for all i <E I , then 
we must have S = {0} or S — IC. 

Proposition III. 4. Let W be a quantum walk according to equations ([I]), ^ and If the coin operators U u on 
which the measure v is supported are such that the set {U*,U U : uj,u' G f2} is irreducible on IC, then W is strictly 
contractive on {I} 1 -, that is, ||W(A)|| < \\A\\ for all A LI. 

Proof. Suppose the set {U^U^ : w, u/ € Vt] is irreducible. It follows from the singular value decomposition of U that 
there exist normalized vectors <f), tp € K. such that 

\\U\\op= {<t>\Uip) = J v(du)(<f>\U u il>) , 

and hence, if ||C/|| op = 1 we must have (filU^) = 1 for all to. This implies (0|C/*/L^0) = 1, hence, 4> is an eigenvector 
of U^Uu for arbitrary uj and ui' , which is forbidden since the set {U*,^ : uj,lu' g f2} is assumed to be irreducible. 



Consequently, \\U\\ op < 1 and by ([21} \\U*.U\\ op < 1. 

Now, let A be {^-independent, i.e. Aq — A, such that A A. 1. Assume ||T(A)| = ||A||, where T denotes the diagonal 
part of C, that is, T(A) = j i/(dui)U*AU u . This implies 

v(du)v(<kj')ti K UZ,A : l\,u:.Al' J .. : tix-.lM. 

On the other hand, we can estimate 

v(du})v(dLj')tT K U*,A*U u ,U*AU u < I I v{du)v{dw')\tr,cU*,A*U u ,UZAU u 



< 



v(du)v{drf)\\XKA\UJJZ,yAU u UZ,\ 
J J v{du)v{<kj')yJ\xK.A* AyJtT K A* A = tv K A* A . 
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It follows from ||T(A)|j = ||yl|| that this inequality is actually an equality, hence, we must have 

\tv K A*(U u U:,)*AU u ,U*,\=tr K A*A, 

which, by the Cauchy-Schwarz inequality, means that A and (U U U*,)*AU U U*, must be proportional for arbitrary uj 
and J . Thus, (U U U*,)*AU U U*, = c^^, ■ A with |c w ^| = 1 and 

The assumption ||T(A)|| = ||^4|| entails c^.^i = 1 for all oj,oj'. Thus, A commutes with U^U*, and since those are 
irreducible it follows that A = a ■ 1. This contradicts the assumption A _L 1, hence, ||T(A)|| < \\A\\. □ 



C. Asymptotic Position Distribution via First and Second Order Perturbation Theory 

Before we start to analyze the first and second order perturbation theory of W e to determine the asymptotic 
position distribution of quantum walks we have to explain why 

W*(l) ->■ nil 

as t -» oo and e -> 0. First of all, we can write 1 = A e + (1 - A e ) and obtain W*(l) = W* (A E ) + W* (1 - A e ). With 
the operator norm estimate ||W e || op < 1, which is valid for e £ K, we get 

l|w*(i)-^ £ ||<|H-^|| 

and consequently W*(l) — > ^* A e with s — ?► and the assertion follows from A e — > 1 as e — >• 0. 

The ballistic respectively diffusive scaling of the position distribution in the asymptotic limit is determined by the 
first respectively second order of the perturbed eigenvalue fx £ = 1 + £// + e 2 /2/i" + . . ., see also }AV WWTT] . And 
indeed, in ballistic scaling, i.e. e = 1/t we get the asymptotic limit of the characteristic function 

C(A) = t Hm M i /t = t lim(l + ^ + ...) t = e-'. (22) 

If, however, the first order is zero // = we may consider the diffusive scaling e = \j\ft of the position distribution 
and get 

C ^ = "Wi = & i 1 + Yt + • ' •) 1 - ^ ■ (23) 

Remark III. 5. In contrast to the case considered in [AVWWlll, where for some models a dependence on the initial 
state po was observed, we see that for the present model the initial state po is irrelevant. 

Now we have all necessary tools to compute the asymptotic position distribution. The first order correction p! can 
be determined from the first order relation obtained from equating coefficients in W e (^4 e ) = p e A £ , which reads 

W'(l) + W(A') = p't + A' , (24) 

where W denotes the derivative of W e at s = 0. The solution to the eigenvector problem W e (A e ) — p £ A e is 
in general not unique. A common choice for A £ is to fix the scalar product of the perturbed eigenvector with the 
unperturbed eigenvector in the following way 

(1|A £ ) = 1 (l|^ n ))=0 ,Vn £ N, 

where A™ denotes the n-th order correction of the eigenvalue A £ . That is, A £ can be expressed as 

n=0 

Fixing the scalar product in this way does not harm analyticity of the A £ , at least for small enough e, since the scalar 
product of an unnormalized vector with 1 is an analytic function in e which is non-zero for small e. 
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The standard approach to determine the higher order corrections to the unperturbed eigenvalue 1 is to expand the 
equation W e (A e ) = fJ, £ A e into a power series in e and then take the scalar product with the unperturbed eigenvector 
1. The choice (t\A e ) = 1, which is equivalent to (l|A( n )) = 0, implies (1\W(A^)) = (1|C(^"))) = for all n > 0, 
thus, some terms in the power series expansion of the eigenvector equation already vanish when taking the scalar 
product with 1. 



By (16) and W £ (A)(p) — S(p)*C(A)(p)S(p + eX) we can express the derivative of W e at e — as 

fvi-X ... \ 
W'(A)(p) = W(A)(p)-iA, A= '■• 

V : Kdim/C'A/ 

which leads to the following expression for the first order correction to the eigenvalue 

M'Hl|W'(l)>^/^^tr K W' (1 ) (25) 



dim/C 

i 

where v is the average of all shift vectors u,. The vector v is closely related to the index indS 1 of the unitary shift 
operator S, see [GNVW12] for a detailed discussion of this quantity. By definition we have 

ind5 = ^Ui (26) 

i 

and consequently v = (dim/C) _1 ind S. Clearly, the characteristic function in ballistic scaling is given by C(A) = e v , 
which is the characteristic function of a point mass at position v. This represents a constant drift of the particle in 
the direction of ind S. It should be noted, that the constant drift is not a ballistic behavior due to the quantumness 
of the walk. Rather, v represents the part of the shift operator, which is not conditioned on the coin state. That 
is, in an extreme case one may consider a one-dimensional quantum walk with two-dimensional internal state space 
JC = C 2 , where vi = V2, i.e. where the shift does not depend on the coin states at all. A ballistic spreading due to 
the quantumness of the walk does therefore not exist in this model. See remark [ill. 6| below for how to investigate the 
diffusive order when v ^ 0. 

Remark III. 6. //ind S ^0 we may subtract the constant drift according to p! — iA • v from the position variable and 
consider the asymptotic distribution of Q = Q — v ■ t. The characteristic function is given by 

C M (A) = trpoW*^^) = trp e- ie * A - 5 W t (e ieA - Q ) , 

and hence we have to consider the modified operator W e (j4) = e~ leX ' v W(Ae 1£X 'Q)e~ leX 'Q = e~ i£X ' v W e (A). The 
eigenvalue fl e is now just given by fi e = e~ leX ' v fi £ , thus 

fie = (l - n'e + ^-e 2 + + f/e + ^-e 2 + 

which shows that ft' = 0. 

In the following we assume f/ = and look at the diffusive scaling of the position distribution. To this end we need 
to determine the second order of the equation 'W e (A £ ) = /J, e A e , which reads 

W"(l) + 2W'(A') + W(A") = n"l + 2fx'A' + A" . (27) 

The second order derivative of the walk operator amounts to W"(A)(p) = W(A)(p) ■ iA = —W(A)(p) ■ A 2 . Again, 
by taking the scalar product with the unperturbed eigenvector we get 

M " = (1|W"(1))+2(1|W'(A')) (28) 

= _j_ r „ i 

(2tt) s 

1 If,. 2i 



/ dSp ^Ar tr ^ (WW + 2iW(A')(p) • A) 
J dim /C 

E "(A • v,) 2 + j^J d* P ^tr, (-W'(l) + A'(p)) • A 



dim/C 
1 

dim/C " %J ' (2ny J ""'dim/C 



' ^(A • v t ) 2 + -L- / 'd s p ^—tr K A'(p) ■ A , 
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where we have used (24) at the fourth step and W(A)(p) = W(A)(p) -iA in every step. Next, we use (24| to eliminate 
A' in this equation. Using // — and W'(l)(p) = iA we get 

W(A')(p) - A'(p) = -iA 

and from p! = 0, which is equivalent to tr^A = 0, it follows that (1|A) = 0. Since the eigenvalue 1 of W is non- 
degenerate we can invert W — id on {I} 1 " and apply it to A e {l} -1 - Denoting the pseudo-inverse by (W — id) -1 we 
get the following expression for A' . 

A'{p) = -i(W- id)- 1 (A)(p). (29) 
Hence, the second order correction pi' can be written as 

«" = ") ! + «bi(hyj d ' pt ' K{w - idr ' wlp> ■ A ' (30) 

i 

Note that the dependence of A on the variable A makes p," a quadratic form in A which defines a symmetric matrix 
D via p" = —X T ■ D ■ A. It is easy to see that the position distribution corresponding to the characteristic function 

C(A) = is a Gaussian with covariance matrix D. 

The following theorem summarizes the results of this section. 

Theorem III. 7. Let W be a quantum walk as defined in ([!]), ^ and ([3|. Suppose for some n e N the power W™ 
is strictly contractive on {I}- 1 . Then: 

i) ballistic scaling: the random variable Q/t, converges to a point measure at v — g^jg indS*. 

ii) diffusive scaling: assuming -5 = the random variable Q/yi, converges to a Gaussian with covariance matrix 
D. The defining equation for D is p" = — X T ■ D ■ A with fj," from (30). Explicitly, the matrix elements of D are 
given by the formula 

2 o f 1 

Daf3 = dhnic^ VhaVh0 + J dSp ^K, tVKRali( " p) ' a '^ = 1 '---' s ' ( 31 ) 

i 

where denotes the a component of the vector Vi and R a p(p) is defined via 

R«p(p) = \(( w - idTHAaXp) • + (W - id)- 1 (A /3 )(p) • A Q ) . (32) 

The diagonal matrices A Q are given by (A Q )y = Vi iCt . 

Of course, for the asymptotic position distribution in diffusive scaling to be well-defined it is necessary that the 
covariance matrix D is positive. This is indeed the case, which we prove now. 



Proposition III. 8. Given the assumptions of theorem III. ? the covariance matrix D is positive, i.e., —fj, > for 
all A. 

Proof. The proof is similar to, though simpler than, the proof of the positivity of the covariance matrix for the models 
considered in [AVWWll . According to (28) we have 

-H" = -<1|W"(1)) - 2(1|W'(A')> 
= (A|A) + 2<iA|W(A')) 

= (A 1 - W{A')\A' - W(A')) + 2(A' - W(A')\W(A')) 
= (A'\A') - (W(A>)\W(A')) , 

where we have used W'(A)(p) = W(A)(p) ■ iA, W"(A)(p) = -W(A)(p) ■ A 2 , iA = A'(p) - W(A')(p) and the fact 
that A' and W(A') are skew-hermitian, which follows from ([29]) and the fact that W(A)* = W(A*) for arbitrary A. 



According to Proposition III. 2 we have ||W|| op < 1, which proves —ft" > 0. □ 



For practical applications it is usually sufficient to give a good approximation to the covariance matrix D. Such an 
approximation can be obtained by a power series expansion of the pseudo-inverse (W — id) -1 . Indeed, if there exist 
n G N such that ||W n || op < 1 we have the following convergent power series expansion for the pseudo inverse. 

(W - id)" 1 = - Wk on the subspace {t} ± (33) 

keN 

This gives us the following corollary, which can be used to approximate the numbers R a p and the covariance matrix 
D. 
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Corollary III. 9. Given the assumptions of theorem III. 7 and v = we have the convergent power series expression 

RMp) = -\Y, ( wfe ( A «)(p) • A £ + W fe (A^)(p) • A Q ) . (34) 

fe£N 

Proof. Let n € N be the smallest natural number such that W" is strictly contractive on {I}" 1 . The convergence of 
the Neumann series for (W — id) -1 can be seen from 

n— 1 oo 

The assertion follows from (1|A) = v = 0. □ 



IV. EXAMPLES 



In this section we apply our theory to five models of randomness and calculate their ballistic and diffusive scaling. 
Continuous as well as discrete distributions on the coin operations are considered and we also include an example in 
two dimensions. In all cases the irreducibility condition of proposition |III.4| is violated, so we have to look at higher 
powers of the walk operator W to show simplicity of the eigenvalue 1. 



A. Coins with Zero Mean U 



In this subsection we assume that v is such that U — J v{du>)U^ = 0. Thus C(A) = C(^4o) = / v{duS)U^ J AoU un 
where we again have abbreviated (2tt)~ s J dp A(p) — A$. If the operator W™ is strictly contractive for some n G N, 



theorem III. 7 is applicable and we can look at the second order correction \i" in order to determine the diffusive 
scaling of the position distribution. Without loss of generality we assume // = 0, see remark [ill. 6[ such that the first 
order of the eigenvalue problem reads W(A')(p) — A'(p) = — iA. Using ~W(A)(p) — S* (p)C(Ao) S (p) , we get 

S*(p)C(A' )S( P )-A'(p)^-iA 

and since A commutes with S(p), this is equivalent to 

CK) +iA = S(p) A' (p)S*(p). 

Now, since the left-hand-side of this equation is independent of p, so is the right-hand-side. Hence, the matrix elements 
of A'(p) are given by (i\A'(p)\j) = aije 1 ^ v '~ Vi ' p with ay € C. Let A' denote the matrix with entries Oy, in other 
words, A'(p) = S*(p)A'S(p). Then, A' can be determined from the p-independcnt equation 

C(P{A'))-A' = -iA, 

where P(A') is defined as 



<W)li> = 



lij , if Vi = Vj 
, if Vi ^ Vj 



B. Hadamard Walk with Random Reflections at Lattice Sites 



A simple example to demonstrate our techniques is the Hadamard walk on t 2 (7L) g ) C 2 which is distorted by random 
reflections. Similar models have previously been studied numerically in |RSA + 05l IPRllj . The intuitive picture is 
that at each time step some links between the lattice points are broken such that the walker cannot pass them but 
is reflected. These reflections can be thought of as flips of the internal degree of freedom followed by the usual shift. 
Hence, the reflections can be implemented by the Pauli operator o~ x which is applied with probability 1 — w for 
W G (0, 1) at each lattice site. 
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The basic Hadamard walk is given by the unitary Wh = S ■ C, where the coin C = 1 ® H and shift operator S are 
defined with respect to a basis |±) of C 2 by 



H = 



1 



y/2 \ 1 -1 



, S\x ® ±) = \x ± 1 ® ±) . 



(35) 



The Pauli operator Uj. acts like a bit flip on |±), i.e. cr^i) = |=p). Following the above, the measure v w on U{2) is a 
convex combination of the point measures Ah respectively A CTx supported on H respectively o~ x . 



v w = w A H + (l — w) A 



(36) 



According to equation ^ the walk operator is given by W(A) = S*C(A)S with shift operator S as in (35) and coin 
operator 



C(A) = ^2 \ x )(y\ ® («> HA XV H + (1 - to) (TzA^o-x) + (1 - 5 xy )UA xy U 



(37) 



■T-.V 



where U = w H + (1 — w) a x . To apply our theory to this particular model, the model has to satisfy the conditions 
of theorem III. 2 Since the set {U*,U U } — {Ha x , a x H} is not irreducible we have to look at higher powers of W in 



order to ensure the non-degeneracy of the eigenvalue 1. 



Lemma IV. 1. Let w G (0,1) and W = S*CS, with C according to (37| and S according to (35|. Then W 2 is 
strictly contractive on {I}" 1 ". 

Proof. We decompose the coin operator again into diagonal part T and conjugation with U, i.e., 

C(A) =T(A ) + U*(A-A )U. 

The operator U = w H + (1 — w) a x is a convex combination of two unitary irreducible operators. Since U is also 
hermitian we get an estimate for its operator norm by considering its eigenvalues. It is easy to see that ||£/|| op = 1 
implies the existence of a common eigenvector of H and o~ x , which is a contradiction to their irreducibility. Hence, 
||l/||op < 1 for all w G (0,1). 

The diagonal part T satisfies the eigenvalue equation T(a y ) = —o- y , hence, ||T|| op = 1. Similar to the proof 
of proposition III. 2 we define the orthogonal subspaces 7o and To according to ( 20 ) . Note that C(7o) C To and 
C(V") C To 1 . It is easy to verify the relations ||W 2 (A)|| = ||C(W(A))|| and \\C(A)f = ||C(A )|| 2 + \\C(A - A )\\ 2 , 
hence, the existence of A _L 1 with ||W 2 (A)|| = \\A\\ implies A,*W(A) £ To- Writing A as a linear combination of the 
Pauli operators A — a x a x + a v a y + a z a z we get the relation 



W(A) = w (a x a z — a y S*a y S + a z S*a x S) + (1 — w)(a x S*a x S — a y S*a y S — a z a z 



(38) 



Hence, ~W(A) € To requires a y = and w a z + (1 — w) a x = 0. Consequently, we have W(A) = (w a x — (1 — w) a z )a z 
and therefore ||W(^1)|| 2 = \wa x - (1 - w) a z \ 2 < {w 2 + (1 -w) 2 ) ■ (\a x \ 2 + \a z \ 2 ) < \a x \ 2 + \a z \ 2 = \\A\\ 2 , which proves 
the assertion. □ 



According to ( 25 1 the ballistic order vanishes, i.e. // = 0. The diffusive order can then be derived from equation ( 30 ) 
and corollary |III.9| The result is depicted in figure[l]and shows that for the completely localized walk, i.e. w = 0, the 
position probability becomes a point measure at the origin whereas for the usual Hadamard walk the diffusion constant 
D diverges. The red points in the figure correspond to an approximation to the diffusion constant by calculating the 
variance of the position probability distributions after 100 time steps. For w < 0.6 the variance data points are in 
good agreement with the diffusion constant, for larger w, however, there is a clearly visible deviation of finite time 
variance and diffusion constant. This is the fingerprint of a generic behavior: when approaching a ballistic quantum 
walk, here w — > 1, the time scale, at which the crossover from ballistic to diffusive behavior happens, diverges. Hence, 
approximations with a fixed number of time steps become worse and worse in the limit w — > 1. 

The plot in figure [l] also shows a comparison of the diffusion constant and the guess w/ (1 — w) by Romanelli et al. 
for the model considered in [RSA + 05 . The diffusion constant for our example shows a similar functional dependence 
on w, though there are small deviations. 



14 




0.4 0.6 

probability w 



1.0 



Figure 1. The diffusion constant of the one-dimensional quantum walk with random reflections as a function of the probability 
w. The blue curve shows an approximation of the diffusion constant evaluated according to corollary |HI.9| for different values of 
w. The red dots show the variances of the position di stribution for t — 100 and the green line shows the guess D — w/(l — w) 
by Romanelli et al. for a similar model considered in |RSA + 05] . 



C. Hadamard Walk with Dephasing 



The Hadamard walk with dephasing is given by a modification of the local coin operator H in Wh — S ■ 1 <g> H , 
with S and H according to ( [35] ) . We modify the coin H by an additional random relative phase shift between the 
states |+) and |— ) 



Note that, by introducing the Pauli operators cr x , a y and a z , we can write H v = e lv<T *H. This phase shift is assumed 
to happen independently and identically in time and space, in other words, the phase ip is chosen for each lattice site 
and in each time step according to a fixed probability measure is on [— tt, 7r). In Fourier space we have the following 
expression for the Hadamard walk with dephasing: 

W(A)(p) = S*(p)C(A)S(p) , C(A) = J v(d(p)He-'^ a ' ■ A ■ e^^H + H*(A~ A )H (39) 

Here, H denotes the operator H = J is (dip) e lv " 7z H . In the following we set r n = | J is(dip)e lvn \ and 9 n = 
arg J is(dip)e lipn and note that C only depends on r n and 6 n for n = 1,2 and not on the actual form of the dis- 
tribution is (dip), see (40) and (41). 
It is easy to see that 

H * iH =H j(<p-?')<T* H = e K v - v ')v x = ( cos(p-p') isin(ip - ip')\ ^ 
v \ isin(ip — ip 1 ) cos(tp — ip') I 



which shows that the set of matrices H^,H V is reducible. Nonetheless, we can apply theorem III. 7 after we prove that 
W 2 is strictly contractive on {l} -1 . 

Lemma IV. 2. Let W be according to ( |39[ ) and v such that r\,r2 < 1. Then W 2 is strictly contractive on {I}" 1 "- 

Proof. We write the coin operator as C(A) = T(Aq) + H*(A — Aq)H and express the dephasing in the maps T 
and conjugation by H in the basis {l,a x ,a y ,a z } of Pauli operators. The dephasing in T corresponds to the matrix 
representation 



A h> / is(d<p)e- iva * ■ A ■ e iipa * 



( 1 o o o\ 

r 2 cos 9 2 —r 2 sin 8 2 

r 2 sin 9 2 r 2 cos 9 2 

\ 1 / 



(40) 
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Figure 2. (a) shows a plot of the diffusion constant D depending on the width 8 of the interval Is from which the phase ip is 
chosen with uniform probability. For the plot shown in (b) we fixed S = n/S and computed the position distribution after 10 
(green), 30 (blue), and 80 (black) time steps with initial state po — 2~ 1 1. For better comparison with the asymptotic position 
distribution (red-dashed line) the average over two neighboring lattice sites was computed such that the finite time step position 
distributions are non-zero everywhere. 



and the dephasing in the conjugation by H has matrix representation 
.1 - / vidipy-**** ■ A- f i/(<V)e i<Aa — > r\- 



( I o o o\ 

cos 26»i -sin 26 1 

sin 29i cos 26 1 

\ 1 / 



(41) 



Clearly, if r\ < 1 the conjugation by H is strictly contractive. If < 1 we have that T is strictly contractive on the 
subspace spanned by a x and a y . The operator a z is mapped to a x by T, and hence, S*T(a z )S G T^, which proves 
||W 2 (A)|| < ||A|| for all A± 1. □ 

Since indS* = 0, the ballistic scaling yields a point measure at the origin as asymptotic position distribution. For 
the asymptotic position distribution in diffusive scaling we need to solve the equation 



W(A')(p)-A'(p) 



which we do by numeric approximation according to corollary |III.9| For concreteness let us consider a family of 
measures vg on [— tt,tt), indexed by 5 £ [0, 7r), and defined via 

Here, Xi s denotes the characteristic function of the set Is = [—5, 6) and dip the Lebesgue measure. The measures us 
represent uniform distributions of (p on Is and figure [2] shows a plot of the diffusion constant and a comparison of the 
position distribution for a finite number of time steps with the asymptotic position distribution. 




D. Quantum Walk with Continuous Coin Distribution 



One of the goals of this paper is to take into account fluctuations of the coin operators in space and time as they 
occur in experiments. However, in most examples up to now, discrete distributions of the coin operators have been 
studied. For a model of realistic noise sources in experiments it seems to be more appropriate to assume a continuous 
distribution around some specific coin operation. Let us therefore consider a quantum walk on ^ 2 (Z) <g) C 2 with the 



shift according to ( 35 ) and coin operator constructed from the following 27r-periodic family of unitary matrices 



U r =( C0Sr Smr , (42) 



sin r — cos r . 
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for which U^/i — H, is the Hadamard coin (35), whereas 11^/2 — &x and Uq — a z . Taking on the idea of random 
fluctuations around some target coin operator, we consider a gaussian probability distribution on /(ro) = [ro — it, 7'o+7r] 
around some point r 6 R. The measure on the parametrization space is then given by 



v(X) = — 



dr, 



(43) 



A" 



where X is a Borel set on /(ro), dr is the usual Lebesgue measure and the normalization N is chosen such that the 
measure is normalized on I(ro). Now we can write the coin operator as 



/ 



C(A)=5»<2,| 



Jx,y 



V I(ro) 



-dre 



\ 



-U*A X yU r + (l-5 x , y )U*A xy U 



(44) 



where 



U 



v{dr)U r — 

I(ro) I(r 



1 (r-rgf 

— dre °- a U r . 



A straightforward calculation shows [U r , U r i] = — 2isin(r — r')a y , hence {U r } is irreducible whereas the set {U*,U r } 
is reducible. Nonetheless, W 2 is strictly contractive on {l} -1 , as the following lemma shows. 



Lemma IV. 3. Let W = S*CS be a quantum walk with C according to (44) and S according to (35), then W 2 is 
strictly contractive on {I}" 1 ". 

Proof. Since the set {U r } is irreducible and all U r are hermitian we can exclude that U has eigenvalues of modulus 
one and it follows that ||E/|| p < 1- This proves the contractivity of the map corresponding to conjugation by U. 
Again, using Pauli operators as a basis for the set of two-dimensional matrices it is easy to see that 



U* aJJ r = — cos 2r <t t 



sin 2r a. 



U r OyU r = 0~y 



U* <7, U r = sin 2r a x + cos 2r a z 



Taking the expectation value of these equations with respect to the probability distribution v shows that the diagonal 
part of C, which we denote again by T, is strictly contractive on a x and a z . The assertion follows from the fact that 
the conjugation with the shift operator S maps a y to Tq~. □ 



The shift operator satisfies indS" = 0, such that according to (25) /jf = which means that the ballistic order 
vanishes. Hence, the characteristic function in diffusive scaling is well-defined and yields a Gaussian. The position 
probability distribution can be calculated by the inverse Fourier transform resulting again in a Gaussian. For a finite 
number of t — 200 time steps the position probability distribution for a coin operator peaked at the Hadamard coin 
= H is depicted in figure |3| (b) for a G {0.01,0.1,0.2,1,2}. It is apparent that for small a the distribution 
looks like the one of the usual Hadamard coin plus some gaussian background. This is again an indication that for 
decreasing c, i.e. coin distributions with decreasing width, the crossover from ballistic to diffusive behavior happens 
at later times. Indeed, since the gaussian measure (43) converges weakly to a point measure, i.e. v — S ro for a — > 0, 



we expect a divergent diffusion constant as the translation invariant and unitary quantum walk with coin U ro exhibits 
ballistic behavior. 

Figure [3] (a) depicts the diffusion constant D as a function of a e (0,1] and r € [0,2n]. Apparently, D < 00 for 
a > independent of the coin at which the distribution is peaked. The symmetry around ro = 7r comes from the fact 
that we have U v+r = a z ll n - r a z . By S z we denote the operator on £ 2 (Z S ) ®JC corresponding to the Fourier transform 
o z . Since Y> z commutes with the shift S we have the^elation J S*E z C(S 2 e ieA Q S z )I] 2 S' = T, z S*C(e ieX ' Q )SE z . Denoting 
the walk operator with additional E z rotations by W, i.e. W(A) = S*Y, z C(Y lz AE z )'E z S we get the relation 



W*(e 



S,W 4 (e UA ' Q )S 



Changing the parameter from r = tt + r to r = n — r is equivalent to changing the initial state from p to S z p ^z- 
Since the characteristic function in the asymptotic limit does not depend on the initial state ( ( [22] ) and (23)), the 
diffusion constant is independent of the above change of ro. The periodic dependence on ro reflects the periodicity of 
the U r . The minima on the slices with a fixed are the ones where the peak of the gaussian measure |43"| ) is at U ro = a x , 
the reason being that the unitary and translation invariant quantum walk with coin a x shows no propagation at all. 
The maxima correspond to ro £ {0, tt} where the measure is peaked around a z . 
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Figure 3. Plot (a) shows the diffusion constant D as a function of a G (0,1] and ro G [0, 2ir]. For every fixed a the diffusion 
constant has a peak at ro G {0, n} (yellow line) which correspond to measures peaked around ±<r z coins. The minima occur at 
ro = 7r/2, Stt/2 (blue line) where the random coins are peaked around a x . The diffusion constant of a quantum walk peaked at 
the Hadamard coin can be found at ro = n/4 (red line) and dephased variants of it at ro = &7r/4, k G 7L. Plot (b) shows the 
position probability distribution of the quantum walk with gaussian distribution peaked at the Hadamard coin H for a number 
of t = 200 time steps in diffusive scaling. Again, the average over two neighboring lattice sites was computed in order to remove 
the oscillating behavior of the probability distributions. The different plots correspond to a = 2 (orange), a = 1 (blue), a = 0.2 
(red), a — 0.1 (magenta), and a — 0.01 (green). What can be seen is that for a — ¥ 0, where the gaussian measure converges 
weakly to the point measure at the Hadamard coin, the diffusion constant diverges, which explains the Hadamard like position 
probability for a = 0.01. Nevertheless for t — > oo it becomes gaussian, as long as a > 0. 

E. Quantum Walks in Two Dimensions 

For a quantum walk on a two-dimensional lattice, the Hilbert space is £ 2 (Z 2 ) ® /C. That is, the position is given 
by two-component vectors x — (xx,X2)- We consider a four-dimensional coin space JC = C 4 = C 2 ® C 2 . The shift 
operator conditioned on the internal degree of freedom is given by 



Sip) = 



( e iv rP 
e lv ^ p 

V 



\ 



P = (Pl,P2) , 



(45) 



where u+ = (1,0), = (—1,0), = (0, —1) and v^. = (0, 1). The explicit example we are going to consider is a 
quantum walk with coin operation constructed from the set of unitaries {Ui = H ® H, U2 = <J X ® 1}, where XJ\ is 
applied with probability (1 — w) and U2 occurs with probability w. In order to apply theorem III.7| some power of 
the walk operator W must be strictly contractive on the subspace {1 <g> I}- 1 . Apparently, W itself is not strictly 
contractive on {1 <g) 1} because the p-independent operator 1 ®i? is an eigenvector of C with eigenvalue 1. Moreover, 
the coins U\ and U2 are reducible. Denoting the eigenvectors of H by |±) it is easy to see that t he tw o dimensional 
subspaces C 2 £g> |±) are invariant subspaces for U\ and L7 2 . To confirm the applicability of theorem 
that W 2 is strictly contractive. 



III.7 



we now prove 



Lemma IV.4. Let w e (0, 1) and W 

and U2 = <J X ® H of which Ui is applied with probability 
strictly contractive on {1 ® 1} . 



S*CS, with S according to (45) and C be defined by the coins U\ = H (g> H 

Then W 2 is 



w and U2 is applied with probability w. 



Proof. Although the coins U\ and U2 are reducible we can exclude the existence of a common eigenvector by calculating 
the determinant of their commutator. If there is a common eigenvector this determinant equals zero, but since we 
have det[J7i, C/ 2 ] = —1 ■ det[_ff, o~ x ] = —2 such an eigenvector cannot exist. This already proves that the operator norm 
of the hermitian matrix U = (1 — w) U\ + w U2 is strictly less than 1, hence, conjugation by U is a strictly contractive 
map. 

Again, let T denote the diagonal part of the coin operator C, i.e., 



T(A ) = (l-w)H ® H ■ A Q ■ H ® H 



W (T, 



1 • A Q ■ <r x 
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By choosing an orthonormal operator basis for To and calculating the action of T with respect to this basis one 
confirms that T is a hermitian map. Clearly, ||T|| op = 1 and in order to verify the contractivity of W 2 we have to 
consider the eigenvectors of T which are orthogonal to 1 ® 1 and have eigenvalues ±1. It is sufficient to prove that 
conjugation by S maps these eigenvectors to vectors with non-zero overlap with the subspace Tq~ ■ These eigenvectors 
have to be common eigenvector of the maps corresponding to conjugation with U\ and U 2 , hence, we find exactly 
1® H, <jy (g> 1 and a y (g> H as eigenvectors with eigenvalues ±1. It is easy to see that the only vectors in 7o which are 
mapped to 7o again are all 4-dimensional diagonal matrices, which proves the assertion. □ 

The argument of the characteristic function is a two-dimensional vector A = (Ai, \ 2 ) T , hence, 



A = 



-Ai 

V 



\ 



A 2 y 



(46) 



Since ^ i v i = 0, the ballistic order is zero by (25), the second order correction p!' can be written as quadratic form 
in A 



-A T • D • A, 



(47) 



with the covariance matrix D. By theorem III. 7 we can compute the matrix elements of D via (31 1 and (32). We use 
again corollary |III.9| to get a numerical approximation of D. 

The asymptotic probability distribution can be computed via the characteristic function 



C{\) = e-^ Tm 

and the probability at position x — (x\, x 2 ) is given by its inverse Fourier transform 

1 



P(x) 



(27T) 2 



K' 2 



dAC(A)e 



— zAa 



(48) 



(49) 



1 



(2tt) 2 



dA, 



1 



2?rVdct D 



see e.g. (NQ98) . with 



D~ 



1 [ D 22 -D 12 
detD \-D 12 D n , 



The asymptotic probability distribution for our example is illustrated in figure [4] for two different values of w. Here 
Ui = H <Ei H renders a Hadamard-walk in two dimensions and U 2 = a x ® 1 yields propagation strictly along the 
diagonal. The coin U 2 in combination with U\ leads to significantly increased spreading in the diagonal direction and 
reduced spreading in the anti-diagonal direction. 



V. GENERALIZATIONS 



Although the vast majority of literature considers quantum walks which are a composition of a single coin and 
a single shift operator we will briefly comment on more general models of quantum walks. A way to characterize 
quantum walks more abstractly is to define them to be discrete time evolutions on a lattice T~L — £ 2 (Z S ) (g> K,, which 
are local and translation invariant. This definition is clearly satisfied if several quantum walks Wj are concatenated 
and considered as a single time step W given by 

W = W„ o . . . o Wi . (50) 

More precisely, we assume that each Wj is a quantum walk according to (pi), which means it can be written as 



W i (A)=S*C i (A)S i , 
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Figure 4. Contour plots of the asymptotic position probability distributions of the quantum walk with coins U\ = H <g) H 
and U<2 — o x ® 1 with probabilities {(1 — w),w} for (a): w — 0.1 and (b): w — 0.9. The coin Ui alone would lead to a 
two-dimensional Hadamard walk, but U2 incorporates propagation strictly along the diagonal. 



where Cj is a, possibly decoherent, coin operator and Si is a unitary state dependent shift operator. For this 
generalized model of quantum walks we have the following proposition, which also covers the extremal case where 
several unitary quantum walks are concatenated with one decoherent quantum walk. 



Proposition V.l. Let W be a generalized quantum walk according to (50). // at least one is strictly contractive 
on {I}' 1 we can apply perturbation theory to the eigenvector equation *W £ (A e ) = [i e A e to determine the asymptotic 
position distribution o/W. In particular, let indSi denote the index of Si according to (26), the asymptotic position 
distribution o/W in ballistic scaling is given by a point measure at (dim/C) -1 X^"=i m dSi- 

Proof. Clearly, W satisfies W(l) = 1 and since W maps 1 and {1} ± to orthogonal subspaces it follows that W is 
strictly contractive on {1}^ if at least one of the W.; is strictly contractive on {l} -1 . Hence, the non-degeneracy of 
the eigenvalue 1 of W is assured. 

The proof of the analyticity of W e is similar to the case in theorem |III.7| hence, W e satisfies the requirements 
of the Kato-Rellich theorem |III.1| This implies that the asymptotic behavior of W can be determined using our 
perturbation method. 

The first order of the power series expansion of ~W £ (A £ ) = fi e A £ reads 

W; o . . . o Wi(l) + . . . + W» o . . . o Wi(l) = //l + A' - W(A') 
and since W 4 (l) = 1 and W' t (A)(p) = Wi(A)(p) ■ iA l this simplifies to 

n-2 

iA n +i^W n o...o Wn-^An-i-Ofo) = (X'l + A'(p) - W(A')(p) . 
i=l 

The scalar product of this equation with the unperturbed eigenvector 1 yields again 

n . n 



□ 



The diffusive scaling of W can also be determined by equating coefficients of the perturbation expansion. Of 
course, the equations get more involved, but the general structure of the problem is the same. For example, the 
second order correction to the eigenvalue fi e is again a quadratic form in A with constant coefficients and the 
asymptotic distribution in diffusive scaling is just a Gaussian independent of the initial state pg. 

To illustrate this, consider a concatenation two quantum walks W = W2 o Wi. The second order equation reads 

W(A") + TW' 2 o W X {A') + 2W 2 o W[(A') + W£(l) + W 2 o W'/(l) + 2W' 2 o Wi(l) = //'l + 2fi f A' + A" 
and by assuming // = and taking the scalar product with 1 again, we get 

which is a quadratic form in A with no further parameter dependencies. 
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VI. CONCLUSION 

We have shown that quantum walks with spatio-temporal fluctuations of the local coin operator exhibit, under 
rather mild assumptions, diffusive behavior in the long-time limit. Our method provides complete information about 
the asymptotic position distribution of the considered quantum walk and, though the appearing equations may not 
be solvable in a simple manner, an approximation to the solution can always be found by computing a truncation of 
a power series expression for the exact solution. 

The model of quantum walks with spatio-temporal coin fluctuations considered exhibits two generic features. First, 
the asymptotic position distribution in ballistic scaling is always given by a point measure. Secondly, for the asymptotic 
position distribution in diffusive scaling we get a Gaussian which is independent of the initial state. 

One crucial assumption in our model is that the coins at different lattice sites and at different times are distributed 
identically and independently. Correlations of the coins in space or time have been studied in the literature and it was 
found that they can lead to different phenomena like ballistic or sub-ballistic behavior or even localization. However, 
a complete theory, combining correlations in time and space and providing sufficiently general criteria for different 
kinds of asymptotic behavior, is still missing. 

The second assumption is of a more technical nature and concerns the irreducibility of the coin operators. This is 
similar to the case considered in |AVWWllj , where it was also shown that the extremest form of reducibility, namely 
commuting coin operators, leads again to ballistic behavior. Our method can, in principle, also be applied to quantum 
walks with spatio-temporal coin fluctuations and reducible coins by using degenerate perturbation theory. 
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